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SIGNIFICANCE  AND  EXPLANATION 

V 

In  our  previous  work  on  predator -prey  systems  we  have  studied  the 
effect  of  constant-rate  harvesting  of  either  species.  The  same  techniques 
may  be  used  to  analyze  the  effect  of  negative  harvesting  rates,  which 
correspond  to  stocking  the  system  by  constant  rate  addition  of  members. 
This  is  a  common  practice  in  many  situations,  particularly  in  connection 
with  fish  populations  where  food  supply  corresponds  to  prey  and  fish  to 
predators,  and  our  resultB  indicate  some  dangers,  principally  the  danger 
of  wiping  out  the  food  supply  by  stocking  fish  too  rapidly.  It  will  be 
seen  that  it  is  safer  to  increase  the  food  supply.  This  avoids  extinction 
dangers,  tends  to  stabilize  the  system,  and  increases  the  equilibrium  fish 
population.  . 
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l.  iotkooktion 

In  out  previous  with  on  predator-prey  ayitMa  (Braurr  and  S»  udack  <1979a>.  hraurr 
and  Soudack  (1979b)  I  w  haw  studied  the  effect  of  constant-rate  harvesting  of  either 
species.  The  sane  techniques  nay  be  used  to  analyte  the  effect  of  negative  harvesting 
rates,  which  correspond  to  stocking  the  system  by  constant  rate  addition  of  members. 
This  is  a  common  practice  in  uny  situations,  particularly  In  connection  with  fish 
populations  where  food  supply  corres|onds  to  prey  and  fish  to  predators,  and  out 
results  indicate  sons  dangers,  pr  inc  l(>al  ly  the  danger  of  wiping  out  the  food  supply 
by  stocking  fish  too  rapidly.  It  will  be  seen  that  it  is  safer  to  increase  the  food 
supply.  This  awnds  extinction  dangers,  tends  to  stabilise  the  system,  and  In.  reases 
the  equilibrium  fish  population. 
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2.  nuatXIKAK'.K.'. 


Me  consider  the  syatr* 

■*  •  -  r 

(1) 

»’  •  yqis.y)  -  C 

<■  a  audel  for  the  sires  *(tl  of  a  prey  population  and  y(tl  of  a  ptedator  population 
at  tin  t.  The  parameters  -r  and  -G  are  non-negative  constants  representing  ratee 
at  which  prey  and  predators  respectively  are  a  {.led  to  the  systen.  Kr  elite  the  atotking 
ratea  in  thia  fora  with  r  ■  0,  C  •_  0  in  order  to  make  the  notation  confora  to  our 
prevloua  work  on  harvesting  (Brauei  and  Soudack  (I9?9a),  Brauvr  and  Soudack  <19?9b>l. 

In  or  ler  to  deacrlbe  the  predator-prey  nature  of  the  system,  we  assume 
(2)  f^ia.yl  <  0,  q^ia.y)  »  0,  g^Ot.y)  ^0 

for  a  *  0,  y  >  0,  and  that  there  raiata  J  »  0  auch  that 
<»  g(J.O)  e  0 

In  aany  of  the  i  aaaarnly-uaed  models.  the  function  q  depends  on  a  only,  correapondinq 
to  the  bioloqical  property  that  the  predatora  do  not  interfere  with  one  another  in  their 
aearch  for  prey.  In  auch  a  aarxtel,  the  curve  q(a.y)  -  0  becomes  the  vertical  atraiqht 
line  a  •  J,  and  q(J.y)  *  0  for  all  y  £  0. 

Nhile  additional  conditions  are  uaualty  lapoaed  (May  (197J),  Pul met  (1976)1,  they 
are  not  required  for  our  analyala.  Under  the  conditions  (2),  (It,  every  orbit  of  the 
ayatea  (1)  with  inlttal  value  in  the  firat  quadrant  of  the  (a.yl  plane  at  t  •  0 
either  leave*  the  firat  quadrant  in  finite  tiae  or  reauilna  in  a  bounded  subnet  of  the 
firat  quadrant  for  all  (  >  0  IBrauer  <19791 1.  Since  F  0  and  G  ^  0,  no  orbit 
can  leave  the  firat  quadrant  aa  t  increases:  therefore  every  orbit  rename  in  the 
firat  quadrant  and  the  poincare-Bendisaon  theoren  la  applicable.  It  followa  that  every 
orbit  tend*  aa  t  *  -  either  to  a  atable  equilibrium  or  to  a  atable  Unit  cycle.  Note, 
-.owe ve r ,  that  there  .  ould  be  several  atable  equilibria  or  Unit  cyclea  with  different 
doaiains  of  attraction.  A  atable  equilibrium  on  one  of  the  coordinate  axea  correaponda 
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to  estlnctlon  of  one  of  the  i|«>  in.  Thin  is  qualitatively  different  biological  )y  fro* 
a  stable  <-  qui  librium  in  the  Interior  of  the  first  quadrant,  oor  rescinding  to  coeslstence 
of  the  two  »|-ocie»,  even  though  there  uy  be  no  qualitative  tutlxiut)  a]  distinction, 
by  analyslnq  the  local  stability  of  equilibria,  we  shall  obtain  information  about  the 
possible  behaviours  of  orbits,  including  estimates  of  the  leqlon  of  asymptotic  stability. 
He  shall  iqnore  the  structurally  unstable  possibility  of  an  equl 1 ibr l us  which  Is  stable 
but  not  asymptot leal ly  stable  on  biological  grounds .  We  shall  a  1 eo  not  consider  further 
the  possibility  of  airs  than  one  stable  limit  cycle,  or  such  situations  as  a  stable 
equilibria*  surrounded  by  an  unstable  periodic  orbit,  which  Is  surroiaided  by  s  stable 
limit  cycle.  Such  poasibl 1  it les  can  occur,  and  our  methods  can  be  adapted  to  deal  with 

tham,  but  we  avoid  them  in  the  Interest  of  siapllctty  of  classl  f  lest  Ion.  However,  we 

shall  see  that  there  ere  situations  In  which  there  must  be  two  equilibria,  each  with  a 
domain  of  attraction  (see  figure  11). 

Since  f  Is.y)  *  0,  the  equation  f(s,y)  -  0  defines  y  ss  a  single-valued 
function  of  s,  denoted  by  y  •  tta),  which  we  may  assvsse  non- nog  at  1  ve  on  soam 

Interval  a  <  I  <  f,  with  i  >  0  and  It  •,  and  f(K,0)  •  0  if  It  «  •.  If  a  •  0, 

It  follows  from  f^is.y)  <  0  that  f(0,0l  >  0,  and  there  estate  L,  0  «  l  <  •,  such 
that 

(4)  f(O.L)  •  0 

biologically,  the  n\aaber  I.  represents  the  meslei*  predator  density  for  which  the  prey 
population  can  establish  Itself  from  e  smell  Initial  population.  If  i  »  0  and 
fta.O)  •  0,  then  f(0,0)  <  Oi  this  Is  the  case  in  which  the  prey  population  Is  unable 
to  develop  even  In  the  ebeence  of  predatore  If  It  gets  too  small.  The  number  f  is  a 
carrying  capacity  for  the  prey  species,  the  miibi*  equilibrium  population  In  the 
absence  of  predators.  In  order  to  Include  such  esanple*  as  the  lake  eutrophication 
model  of  'brlen  (197«i  and  the  spruce  budvem  model  of  t.udwlg,  Jones,  and  Helling 
(1978),  we  permit  the  possibilities  K  •  •  and  l.  -  • . 
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The  qualitative  behaviour  of  the  system  (1)  la  quite  different  in  the  three  ceaea 


(ll  0*0,  (It)  a  .  o  but  l  •  •,  and  (tit)  a  •  0,  l  *  •.  Accordingly,  It  will 

be  necessary  to  dlatlngulah  theee  caaea  in  what  tollowa. 

Since  g^lx.y)  »  0  and  g^lx.y)  ^  0,  the  equation  g(x,y)  *  0  definea  a  aa  a 
monotone  non-decreasing  function  of  y,  denoted  by  x  •  Hy).  In  view  of  the  assunp- 
tlon  (2),  thia  function  la  defined  for  0  •_  y  *  •,  and  J  •  CIO).  In  aiany  of  the 
■Odets  used  to  describe  predator-prey  systems.  the  function  g  la  independent  rt  y, 
aa  w«  have  already  Mentioned.  In  this  case  T  is  the  constant  function  a  •  J. 
Biologically,  J  la  the  minimi*  prey  population  required  for  the  predator  population 

to  establish  Itself.  Me  will  examine  explicitly  only  the  cese  J  <  K;  the  case  J  ^  I 

can  be  treated  in  exactly  the  same  way,  with  analogous  resulta. 

An  equilibrium  of  (1)  la  an  Intersection  <a,y)  of  the  prey  Isocline 

(5)  af(a.y)  -  r  •  0  , 
and  the  predetor  isocline 

(6)  ygla.y)  -  C  ■  0  . 

The  assimifttona  we  have  made  on  t  and  g  ig  ly  the  existence  of  an  equilibrium 
Ix^.y^)  with  *m  »  0,  >  0  for  f  •  0,  C  •  0.  Me  shall  asaimw-  that  this 

equilibria,  in  the  interior  of  the  first  quadrant,  Is  unlqus.  Aa  r  and/or  C  vary, 
the  equilibria  (x_,y^)  will  move,  and  nay  disappear  for  some  critical  stocking  rates. 
If  slther  r  or  C  Is  rero,  there  may  also  be  equilibria  on  the  coordinate  axes. 

The  slope  of  the  prey  isocline  et  en  equilibrium  (x,y)  is 

xf  (x,yl  ♦  f (x,y 1 
xf  <x,y) 

which  has  tha  same  sign  as  xf^lx.y)  ♦  f(x.y),  by  (2),  while  the  slope  of  the  predator 
isocline  is 

Y^^fx.y) 

igy(x,y)  «  g(x.y) 

which  has  opposite  sign  to  yg^lx.y)  ♦  g(x,y)  by  (2)  unless  g(x.y)  •  0  (that  Is, 

C  •  0)  end  g^lx.y)  •  0,  )n  which  case  the  predator  Isocline  is  vertical.  By 
•umlnh  i  the  coefficient  matrix 


Mx.yl 


Ixf^lx.y)  ♦  f(i,y)  xfy<«.y) 

yq^lx.yl  yqy(i.y)  ♦  q<i,y) 

o(  the  linearised  system  at  an  equilibrium  (x,y)  and  observing  that  <X,y)  is  a 
saddle  point  >f  and  only  it  dot  Mx.yl  <  0,  it  la  eaay  to  verity  that  <x,y)  la  a 
saddle  point  if  and  only  If  the  alope  of  the  predator  laocllne  la  either  negative  hut 
larger  than  the  alope  of  the  prey  laocllne  or  poaltlve  but  aouller  than  the  alope  of 
the  prey  laocllne.  It  the  predator  laocllne  la  vertical,  (ir.y)  can  not  be  a  saddle 
point.  This  generalizes  a  criterion  of  Hulmer  (1916)  for  F  *  0,  C  •  0.  In  particular, 
the  remark  of  Bulaetr  that  If  there  la  a  unique  equilibrium  in  the  interior  of  the  flrat 
quadrant,  then  thla  equilibrium  can  not  be  a  saddle  point,  eatenda  to  arbitrary  F 
and  C. 

In  the  follovinq  aectiona,  we  shall  examine  the  effect  of  varying  F,  or  C,  or 
both  F  and  C  simultaneously.  We  shall  pay  particular  attention  to  the  matheraat leal 
catastrophes  which  may  arise  if  the  equilibria  which  la  at  (x^,y^ I  for  F  •  C  •  0 
coalesces  with  the  equilibrium  which  is  at  the  origin  for  F  *  C  •  0.  There  la  another 


saddle  point  which  is  at 

(K.0) 

for 

r  • 

C  -  0,  but  we  shall  see  that  It  plays  no 

role  in  the  analysis  when 

F  «  0 

or 

c  « 

0. 
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J.  PKEY  STOCKING 


*1«  G  «  0  and  Juitrisc  r  fro*  sero.  If  a  •  0.  then  for  E  •  0  there 
»•  "»  equilibria  Which  l*  nor  a  saddle  point  and  there  are  two  saddle 

points  S^(0,0)  and  S^IK.O).  If  a  »  0,  then  for  t  •  0,  Sj  la  an  asymptotically 
■  table  equilibria,  and  there  la  a  second  saddle  point  E^la.O).  As  T  decreases, 
the  prey  isocline  Moves  up,  since  fy(s.y)  *  0.  Tor  r  <  0  the  prey  Isocline  Is 
asymptotic  to  the  y-asls  and  to  the  curve  f(s,y)  -  0.  Thus  if  a  -  0,  sws 

into  the  second  quadrant  and  S Moves  to  the  rlqht  alonq  the  s-asis  (figure  1),  If 
a  '  0,  Sj  M)V»I  to  the  right  and  *>vea  to  the  left  alonq  the  s-asls  until 

and  5^  coalesce  and  disappear  (Eiqurr  J). 


g(x,y)=0 


J 


- SADDLE  ASYMPTOTE 

U  REGION  OF  INITIAL  STATES 

LEADING  TO  PREDATOR  EXTINCTION 
AT  Si 


1 


s2 


Since  every  orbit  which  at  arts  in  the  first  quadrant  at  t  •  0  remains  In  a 
bounded  subset  of  the  first  quadrant  for  all  t  >  0,  the  Poincare-bendixaon  theorem 
shows  that  every  orbit  tends  either  to  P^  or  to  a  limit  cycle  around  p^,  or  to 
(If  a  >  01  as  t  •  with  the  esceptlon  of  orbits  start inq  on  the  x-axis  and 
tendinq  to  the  saddle  point  and.  in  the  case  a  •  0,  a  saparatrlx  in  the  first 

quadrant  tendinq  to  the  saddle  point  (Figure  2).  In  the  case  a  »  0  this 

separatnx  divides  the  first  quadrant  into  two  regions,  the  dosiain  of  attraction  of 
P^  or  a  limit  cycle  around  P  ,  for  which  both  species  coexist,  and  the  domain  of 
attraction  of  Sj,  for  which  the  predators  become  extinct  as  t  ♦  •. 

As  F  decreases,  the  equilibria  P^  serves  up  along  the  curve  q(x,y)  •  0.  If 
sf^ls.y)  ♦  f(x,y)  «  0  for  larqe  y,  as  is  the  case  for  the  coassonly  used  predator- 
prey  models,  then  It  Is  essy  to  see  thst  the  rqulllbnum  P^  must  be  asymptct  leal  ly 
stable  for  all  larqa  F  ,  This  suggests  that  If  a  •  0  there  are  essentially  two 
possibilities.  Either  Pm  is  aaymptot  1  cal  ly  stable  for  all  T  0  and  all  orbits 
tend  to  P^,  or  P^  Is  unstable  initially  and  orbits  tend  to  a  limit  cycle  around 
P  but  aa  r  increases.  P  stabilises  and  orbits  tend  to  P  . 

If  a  >  0,  the  behaviour  in  the  domain  of  attraction  of  P^  or  a  limit  cycle 
around  P_  Is  the  same  as  for  the  case  a  *  0.  As  F1  Increases,  the  domain  of 
attraction  of  shrinks  until  and  coalesce.  This  coalescence  is  a 

mathematical  catastrophe,  though  without  the  usual  connotation  of  disaster.  It 
signals  s  transition  to  a  qualitative  picture  like  the  rase  for  a  •  0  with  no 

possibility  of  predator  extinction)  there  exists  F  <0  with  predator  extinction  for 

c 

some  Initial  states  If  0  >  F  -  F  and  coexistence  for  all  initial  states  if  F  •  F  . 

—  c  c 
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4.  t-MlVATt'K  STOCKING 


We  mm  fia  r  •  0  UK)  decrease  C  turn  into.  For  C  •  0.  the  |  ■  odator  isocline 
(6)  la  asyuptot  ic  to  the  curve  gla.y)  •  0  aa  y  •  “  and  Ilea  to  the  left  of 
gla.y)  -  0  because  of  (I).  In  order  to  analyse  the  behaviour  aa  G  decreases  we  eust 
>«1  Jiato  the  thiee  posslbl  1 1 1 leal 

la)  n  ■  0  eo  that  the  curve  fla.y)  •  0  does  not  intersect  the  |»eltlve  y-aala 

(bl  i*0  but  L«e,  so  that  f(0,y)  »  0  for  y  »  0  amt  I  la  flat  • 

■  •O, 

Ic)  »  •  0  and  there  eaiata  L  *  *  with  f(0,L)  •  0. 

If  >  '  0,  as  G  decreases  r_  reeves  along  fla.y)  >0  to  the  left.  toves 

ui-  the  y-aais.  and  a  vc i  along  fla.y)  «  0  to  the  right,  while  ra  ves  down 

into  the  f  urth  guadrant  (Figure  1).  Just  as  in  the  |  rey  stocking  aituatlan,  there 
la  a  seiaratna  at  which  dl vitas  the  first  guadrant  Into  a  detain  of  attraction 

of  or  a  Halt  cycle  aroiaid  1_  for  which  both  species  consist,  and  a  dcasaln  of 

attraction  of  S^.  The  only  difference  Is  that  now  S(  oorres|>nnds  to  prey  eattnctlon. 
eventually.  and  t#  coalesce  and  dleapt«ar.  When  this  occurs,  we  have  a 

aat  hesaat  leal  I  ataatrophe,  amt  all  orbits  tend  to  (prey  eatinctlan).  Thus  there 

eatsts  G  «  0  with  consistence  for  eoae  Initial  states  If  0  >  G  »  G  and  prey 

c  —  C 

eattnctlon  for  all  Initial  states  if  G  *  G  . 

c 
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The  situation  in  which  a  -  0  and  L  *  *  la  considerably  tort  cuau  l  Rated. 
Slnca  this  is  the  one  which  describes  an  at  of  the  cuaaonly  used  aixlels,  we  esaaine 
It  in  detail.  As  C  decreases,  aoves  to  the  left  along  f<».yl  •  0  and 

aivti  up  the  y-asis,  with  every  orbit  starting  in  the  interior  of  the  first  quadrant 
tending  either  to  P_  or  to  a  limit  cycle  around  as  t  4  •  (figure  SI. 


he  Sj  reaches  the  Point  Q<0,L)  where  f(B,y)  •  0  seats  the  y-asls,  that*  arc  two 

possible  caaaa: 

Case  1 1  The  slope  of  tha  pray  Isocline  fo>  f  *  0  la  leas  than  tha  slope  of  tha 
predator  laocllna  with  C  chosen  so  that  tha  predator  laocllna  paaaaa  through  10,  L) 
(Pleura  4).  In  thla  cane,  Sj  and  p^  aaat  at  0-  *«fore  thla  occurs,  every  orbit 
tends  to  P  or  to  a  ll»lt  cycle  around  P  as  t  *  ha  -  G  Inciaasas  further, 

Ps  disappears  and  becomes  an  asymptotically  stable  node  to  which  every  orbit 

tends  as  t  •  corresponding  to  prey  eatlnctlon.  Soanhst  as  in  the  situation  in 

which  a  »  0,  there  ealsts  C  with  consistence  for  all  initial  states  if  0  »  G  »  G 

C  "  c 

and  prey  eatlnctlon  for  all  initial  state  If  G  *  Gc-  The  critical  stocking  rate  -Cc 
is  qiven  by 

(7)  -Cc  -  -Uj(O.L) 
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ceee  1 1 The  elope  of  the  prey  let*  line  for  K  •  0  u  greeter  then  the  elope  of 
the  predator  laocllne  with  G  choeen  eo  thet  the  predator  leocllnr  paeeea  through 
(0.  L)  (Figure  7).  In  thle  ceee  S  ^  r  cache  e  (J  while  P_  it  etill  In  the  Interior 
of  the  firet  quadrant.  Until  thle  occure,  every  orbit  tende  to  or  to  e  limit 

cycle  eround  P^  ee  t  •  *.  When  paaaee  Q,  e  new  equilibria  eppeere  -  e 

eeddle  point  T  on  fla.y)  •  0.  The  eteble  eeperetrlcee  et  T  separate  the  firet 
quadrant  into  e  region  for  which  orblte  tend  to  P_  or  e  limit  cycle  eround  P_  end 
a  region  for  which  orblte  tend  to  the  eteble  node  S^(  e  region  of  prey  extinction 
(Figure  (I).  There  ie  e  eecond  treneltlon,  when  T  end  P^  roeleece  end  dleeigieer, 
end  ell  orblte  tend  to  (Figure  9).  Thue  in  thle  ceee  there  ealet  Gc<  G* <  0 

with  coealetenco  for  ell  Initial  atetee  if  0  ^  G  »  Gc>  e  region  of  coeatetence  and 
e  region  of  prey  extinction  if  Gc  >  G  »  C*.  end  |>rey  extinction  for  all  Initial 
etetee  If  G  <  G* .  he  In  Ceee  1,  Gc  la  given  by  (7),  and  la  eaetly  calculated.  The 
eecond  treneltlon  rate  G*  can  aleo  1*  calculated  ee  the  value  of  G  for  which  the 


interaect lone  of  (6)  with  fla.y)  -  0  coaleace. 


L 


Fig.  8 
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Since  the  distinction  between  ase*  I  snd  II  is  In  t«r»  of  the  slopes  of  the 


isoclines  it  (O.L),  It  is  easy  to  see  thet  the  systea  (1)  is  in  Csse  I  if  end  only  If 

f  (O.L)  Lg  (O.L) 

(8)  -  -5 -  <  .  - 5 - 

f  (O.L*  4(0. L)  ♦  Ut  (O.L) 

y  v 


For  the  1 at  ge  class  of  Models  elth 

f («,y)  •  «(s)  •  yh(s) 


(9) 


4(s,y>  •  slsh(s)  -  Jb(J)) 

where  » (s)  0  for  0^*«K,  ♦  '  U)  _  0,  htslj^O,  h’(s)^0,  [xh(s))'>0. 

and  *h(s>  is  bounded  «s  s  •  -  (Maynard  Salth  (1974)),  the  condition  (8)  reduces  to 
Jh(JI|h(0l*' (0)  -  «(0)h'(0>)  •  •<0)(h(0)|]  . 
frees  which  we  see  that  a  sufficient  condition  for  Case  I  Is 


(10) 


Jh(Jl 


(h(0>)2 


*  1 


(Observe  that  h(0)  '0  if  L  «  •).  It  is  easy  to  show  that  If  |  -•  ^  - 1  0,  then 

(10)  Is  satisfied,  and  therefore,  so  are  the  standard  forsts  of  predator  functional 

response  such  sa  h(s)  -  — — — -  (Moiling  (196%))  and  h(s)  •  -  (1  -  e  CB)  llvlev 

x  •  A  X 


(1961)). 


However,  It  Is  not  true  that  all  asodcla  of  the  font  (9)  are  in  Case  t. 
(a-ksittedly  artificial)  choice 


h(s) 


1  -  (1  ♦  s) 


-1/2 


The 


for  which  h(0)  -  h*(0)  •  -  ,  violates  (10)  If  J  >  8.  More  significantly.  In 

Section  6  we  shall  describe  a  class  of  predator-prey  node  I  a  In  which  predators 
interfere  with  one  another  for  which  the  systea  II)  is  In  Case  11  if  this  interference 
is  large  enough.  Thus  the  biological  effect  of  predator  Interference  Is  to  give  s 
possibility  of  prey  survival  (depending  on  the  Initial  atate)  under  predator  stocking 
which  would  guarantee  prey  estinctlon  without  this  Interference. 


If  a  *  0,  the  situation  Is  slightly  different  when  -r  is 


Her  than  the 


stocking  rate  fur  which  ami  coalesce  (Tig ure  ill.  There  are  tsu  nenslrhit 

lefluns  divided  by  the  stable  sspatatrlces  t end l ml  to  the  saddle  Joint  S  One  is 
the  dasin  of  attraction  of  Pm  or  a  limit  cycle  around  p  .  and  the  other  Is  the 
dtssain  of  attraction  of  the  asymj>tottcal  ly  stable  e>|ul  1  Ibr  1  un  8  ,  wtilch  is  in  the 
Interior  of  the  first  quadrant.  The  ooeslstance  teuton  is  the  entire  first  quadrant, 
but  there  are  two  possible  equilibria,  depending  on  the  initial  state. 

There  are  acany  practical  questions  whuh  can  be  raised  for  tw>>-spectes  st  .  king 
problems.  Tor  sras) le,  what  is  the  effect  of  increasing  one  stocking  late  while 
holding  the  other  flaed/  Such  questions  may  be  studies  met  readily  by  a  output  er 
simulation.  In  the  following  section,  we  shall  Indicate  the  lesults  of  sosw  simulations 
for  a  class  of  eias| les. 


L  // 


b.  tXAMnX.-. 


Ill) 


i»a  cmaaonly  uaad  n«l»l  for  prodator-proy  intn^Uuni  I  Hoi  1  wxj  U4bS)|  ta 

fla.yl  -  r(l  -  J)  - 
,(a.yl  -  -*U- 


la  •  A)  IJ  ♦  A) 

It  nay  aaaily  ba  iMiirxl  (toa  our  aarltar  atudlaa  on  lurwatlnq  lluuri  and  Soudack 
(1974a),  Hi  ajar  and  Soudack  (1979b)  1  that  undar  pula  pray  atock  ing ,  tha  aqulllbrlv» 
»»  a»ya|  toll,  rally  a  tab  la  If 


-r  » 


rJ 


IK  -  A  -  2J)  . 

u 

Undar  pura  pradator  atockirvj,  tha  oquillbrlva  la  aaynptot leal ly  atabla  l( 
ra^U  ♦  A)  liMm  ♦  A  -  K)  ♦  aAKIJ  •  |J  >  0  | 

tha  corraaiondlnq  a took ing  rata  nay  ba  caleulatad  (run 

r*A 


IMa 


-®  -  iTTJTTr  M  '  *  •-» 

1  la  of  tha  typa  a  “  0,  L  <  •  ,  with  L  •  rk.  Tha  pradator  atocklnq  rata 


for  which  tha  aaddla  point  raachaa  tha  point  Q(0,L)  ta  qlvan  by 


(13) 


raAJ 
-C  •  — 
C  J  ♦  A 


Aa  wo  hav*  ahown  In  Sac t l on  4,  thla  nodal  la  nacaaaartly  In  Claaa  I. 

A  variant  of  tha  nodal  (11)  which  incorporataa  tntarfaranca  aennq  tha  pradatora 
haa  boar,  auqqaatad  by  Aatykin  (1974).  Por  thla  andal 


.y)  *  r(l  *  f )  *  pf-j 


(111 


f  (a 


gla.y)  • 


aAla  -  J) 


-  wy 


(a  •  A)  (J  •  A I 

ho  paranatar  w  nay  ba  vlawad  aa  a  naaaura  of  tha  anount  of  intarfaranca.  Thla  nr>dal 
la  alao  of  tha  typo  a  •  0,  L  «  •  with  L  •  rA,  but  now 


(14) 


_  raJA  „  3  3 


'oaparltnn  with  113)  ahowa  that  for  thla  andal  a  htqhor  rata  of  pradator  atockinq  la 
naadad  bafora  a  tranattlon  la  raachad  than  whan  thara  la  no  pradator  Intarfaranca. 


By  eaanining  the  condition  (B),  «  observe  that  the  Model  til)  la  in  Caae  1  if  and 


only  If 
(1S) 


to 


a  IK  «  J) 

Jr l*  -  Al  1J  ♦  A) 


Thus  there  la  a  predator  atockirw  range  in  which  either  both  aperies  coealat  or  the  prey 
a  pec  lea  bcaaea  eatlnct.  depending  on  the  Initial  state,  whenever  there  la  enough  inter¬ 
ference  easing  the  predators  In  searching  for  prey. 

This  principle  generalises  ivedlately  to  Models  of  the  for* 
f(a,y)  *  «(a)  -  yhia) 
g(a.y)  •  stahla)  -  Jh(J)|  -  i»o(y)  , 

with  «(a)  2.  0  for  0  L  *  L  *•  ♦'*■>!  0,  him)  ^  0,  h'  (a)  ^  0,  tahlair  »  0. 
ah  (a)  tnuxlad  as  a  *  •,  t(0)  •  0,  o'  (y)  »  0,  s'  (0)  *  1.  ft  Is  easy  to  see  that 
such  s  aside  1  la  In  Case  I  If  and  only  If 

|*‘  (0)h(0)  -  tt0)h*  (0)1 (sJh(J)  ♦  w<a(L>  ♦  a*  (L)  1  *  •♦<0><M0)|J 

nioi 

Since  a(L)  *  a*  (LI  l >  0,  this  condition  can  be  violated  by  taking  to  sufficiently 
H  (01 

large.  Thus  the  Model  111)  Is  of  the  type  a  *  0,  L  «  •,  and  May  be  either  Case  1 
or  Case  II. 

Another  variant  of  (11),  based  on  s  suggestion  of  Clark  lUTk),  p.  1M,  la 
f(a.y)  e  r(a  -  a)  (  1  - 


(16) 


g(a.y) 


-  J  l _ 

(a  •  A)  (J  ♦  A) 


It  is  easy  to  see  that  this  la  of  the  type  a  »  0.  The  spruce  budwora  node)  of  Lwbrig. 
Jones,  and  Holllng  (19J8)  Is  of  the  fora 

f (a.y)  •  rd  -  |l  -  — j 

a  *a 

gla.y)  -  sU  -  ;£) 

which  Is  easily  seen  to  be  of  the  type  a  •  0,  L  •  •.  Thus  all  the  types  of  Model 
we  have  considered  are  realisable  In  practice. 


-J*' 


Ms  have  (silotaiJ  some  computer  simulations  for  th*  l  (11)  with  r  •  1,  •  •  |, 

K  *  40,  J  •  10,  A  -  10,  It  is  sssy  to  calculate  that  the  equlllbritmi  (a  ,y  )  is 
a»y*p  tot  teal  ly  stable  uiilri  pure  prey  stocking  if  and  only  if  -r  A  SO  and  under 

pure  predator  stocking  if  and  only  if  -G  »  1. 1871,  Further,  -C  -  S. 00,  and  the 

c 

variant  (11)  with  the  iaar  parameters  is  in  Case  11  for  ti  »  If  s  *  0.0%,  then 

-C^  -  10.00  for  (11).  These  results  can  of  course  be  contirrsed  by  cosgruter  simulation. 

Computer  simulation  is  s  is  useful  in  esaming  the  behaviour  of  the  system  under 
simultaneous  stocking  of  both  species.  ror  esample,  under  Used  prey  stocking, 

-f  •  1,  we  found  that  for  -O  •  0.5  there  is  a  limit  cycle,  which  heomes  an 
asymptotically  stable  spiral  point  for  *C  •  1,  and  an  asymptotically  stable  node 
for  -C  •  4.  However,  as  -C  increases,  the  equilibria  approaches  the  y-asis, 
raising  the  possibility  of  prey  estlnction  caused  by  a  small  perturbation 
(Figures  12,  11,  14). 
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7. 


MTUISlOMh 


It  appears  that  qualitatively  the  effect  of  increasing  one  of  the  stocking  rates 
while  holding  the  othet  one  fixed  ix  quite  sisular  to  the  effect  under  xinqle-xpeciex 
stocking  .  Predator  stocking  tend*  to  proaote  prey  extinction,  which  could  have 
practical  appli.  at  ion*  fur  a  |>e»t  1  lent  lal  prey  species,  provided  the  predator  specie* 
either  require*  prey  to  survive  or  1*  itself  not  a  nuisance.  Prey  atockinq  tend*  to 
stabilize  the  system  and  increase  the  predator  populat Ion,  which  could  have  practical 
ap|lications  for  a  ptclator  specie*  which  one  wishes  to  harvest  in  quantity. 

If  the  predator  specie*  is  a  fish  and  the  prey  specie*  Is  Its  food,  our  results 
indicate  that  indiscriminate  fish  stocking  may  exhaust  the  food  supply,  leading  to  a 
collapse  of  the  fish  population.  Food  svfply  stocking,  possibly  coupled  with  Modest 
fish  stocking  Is  a  safer  and  anre  productive  procedure.  Onaqmter  emulations  may  be 
of  great  use  in  determining  optimal  stocking  rates,  taking  Into  account  not  only  such 
factors  as  stability,  but  also  food  coats.  Our  results  Indicate  some  directions  to  be 
taken  in  analyzing  mixed  xt.  eking  and  harvesting  situation*,  which  could  have  practical 
significance.  We  propose  to  continue  with  this  investigation. 

The  authors  wish  to  thank  Mr.  A1  Mackenzie  for  drawing  the  nisseroua  figures 
which  apf>ear  In  this  paper . 
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